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With the help of a conformal, timelike Killing-vector we define generalized equilibrium states 
for cosmological fluids with particle production. For massless particles the generalized equilibrium 
conditions require the production rate to vanish and the well known "global" equilibrium of standard 
relativistic thermodynamics is recovered as a limiting case. The equivalence between the creation 
rate for particles with nonzero mass and an effective viscous fluid pressure follows as a consequence 
of the generalized equilibrium properties. The implications of this equivalence for the cosmological 
dynamics are discussed, including the possibility of a power-law inflationary behaviour. For a simple 
gas a microscopic derivation for such kind of equilibrium is given on the basis of relativistic kinetic 
theory. 
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(N| ' I. INTRODUCTION 

Simple relativistic gases or fluids may be in a state of "global" thermodynamic equilibrium if the spacetime admits a 
timelike Killing- vector (KV) i.e., if it is stationary Jl] — p~0|| - The immediate consequence of this well-known fact is that 
qq ' such kind of equilibrium cannot generally exist in the homogeneous and isotropic Friedmann-Lemaitre-Robertson- 
Walker (FLRW) models of the expanding Universe since these models are not stationary. The only exception is a 
7— I ■ system of relativistic, massless particles (radiation), for which the "global" equilibrium condition is less restrictive 
and requires the existence of a conformal timelike Killing-vector only. A conformal timelike Killing-vector however, 
is admitted in FLRW spacetimes. 

As was pointed out recently |]ll| , |l2| , this general picture representing common wisdom is modified if one drops the 
CJ assumption of particle number conservation which is always implied in the mentioned standard considerations. If 
the number of simple gas particles is increasing, i.e., if there is net particle production (of whatsoever origin), there 
Ch . exist collisional equilibrium states which we will call here "generalized" equilibrium states, relying on the existence 
of a conformal Killing- vector (CKV), for any equation of state admitted by the standard kinetic theory of a simple 
, relativistic gas, not only for radiation. Radiation under these conditions is the limiting case of vanishing particle 
production. Denoting the fluid energy density by p, the equilibrium pressure by p, the temperature by T and the 
fluid four- velocity by u l , it was shown that for any equation of state within the range between nonrelativistic matter 
(p <C p) and radiation (p = p/3) there exists a well-defined rate at which the corresponding particles have to be 
produced in order to satisfy the (generalized) equilibrium conditions following from the CKV property of Ui/T. This 
rate is highest for nonrelativistic matter (p -C p) in which case it amounts to half the expansion rate and it vanishes, 
as already mentioned, in the opposite limit of radiation ( p = p/3). 

All results concerning this kind of "generalized " equilibrium states characterized by a timclike CKV under the 
condition of increasing particle number were obtained within a simple particle creation model on the level of relativistic 
kinetic theory. A modified Boltzmann equation was proposed in which a "source term" describes the change of the 
one-particle distribution function due to particle number nonconserving processes, supposedly of quantum origin. The 
above statements on the creation rates necessary for the generalized equilibrium conditions to be fulfilled were then 
obtained from a simple rate approximation of the source term in Boltzmann's equation. Therefore the question arises, 
whether the new equilibrium states associated with a varying particle number are just features of a specific creation 
model or whether they reflect a more general property of corresponding physical systems. 

The main purpose of the present paper is to clarify this situation. It will turn out that the characterization of 
generalized thermodynamic equilibrium through a conformal, timelike KV for a fluid with nonconserved particle 
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number as sketched above is completely general and does not depend at all on a specific creation model. This 
conclusion is obtained on a purely phenomenological basis without any assumptions other than standard ones in fluid 
cosmology, except for admitting a nonvanishing source term in the particle number balance. 

Furthermore, the equivalence between a nonvanishing cosmological particle production rate and an effective bulk 
pressure is shown to be a consequence of the generalized equilibrium conditions provided standard assumptions for 
imperfect fluids are made. Especially, it will not be necessary to introduce a "creation" pressure by hand. This 
supplementary pressure may be derived and needs not to be postulated. 

The second purpose of the paper is to overcome some limitations of the previously used kinetic model for particle 
production QJl^l • Within the latter model the equivalence between production rate and effective bulk pressure was 
restricted to spatially homogeneous spacetimes. We show here that this restriction is unnecessary. We will establish 
an exact correspondence between the kinetic theory model for particle production and the general, model independent 
phenomenological approach. 

The plan of the paper is as follows. In Sec. II we define a phenomenological concept of generalized equilibrium 
which is compatible with an increase in the fluid particle number. Section III investigates the thermodynamics of 
the corresponding creation process with special emphasis on the limiting cases for pure radiation and nonrelativistic 
matter. Section IV considers the cosmological dynamics for a universe filled with matter in generalized equilibrium. 
Section V presents a microscopic justification of the generalized equilibrium concept for a relativistic gas. A brief 
summary of the paper is given in Sec. VI. Units have been chosen so that c = &b = h = 1. 



II. GENERALIZED EQUILIBRIUM 

The energy-momentum tensor of an imperfect fluid is generally given by (see, e.g., |8|) 

pik 
(0) 



Tik rr\ik , j ik , ik , i k , k i /-i \ 

= 1 (0) + irn + 7r + q u + q u (1) 



with 



I$j= pu l u k +ph ik (2) 



and 

ik i i iik n i -i /o\ 

ir Uk — q Ui = 7Tj = h Ui = I) , u Ui = — 1 . (6) 

Here, p is the energy density of a fiducial equilibrium state, p is the corresponding equilibrium pressure, u l is the fluid 
four-velocity in the Eckart frame and h ik = g %k + u l u k is the spatial projection tensor. The quantity ir denotes that 
part of the scalar pressure which is connected with entropy production, — ir lk is the anisotropic stress tensor and q l is 
the heat flux vector. Within the Eckart frame the particle number flow vector N l is given by 

N l = nu l (4) 

where n is the particle number density. While we require local energy-momentum conservation T lk k — in accordance 
with the integrability conditions of Einstein's equations, we will not assume the fluid particle number to be conserved. 
Instead, we admit a source term T in the corresponding balance equation which describes the rate of change of the 
fluid particles: 



Inserting here the expression (H) yields 



m = nT , (5) 



h + On = nT (6) 



where O = u l . { is the fluid expansion and h = n.iU 1 . The quantity T is assumed to describe the net change of the 
particle number either due to internal reactions within the medium ]T^Jl4] | or due to particle production in strongly 
varying gravitational fields H-0]. Usually, it is an input quantity in a phenomenological description which for 
specific processes has to be calculated from the underlying microphysics. In our considerations T will be fixed by the 
generalized equilibrium conditions to be discussed below. 

Projecting the covariant derivative of T lk in direction of Ui provides us with 

u{T! k = Ul Tl k yk + Ui [7Th ik + ir* k + q*u k + q k u>] . fc = . (7) 
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The first term on the right-hand side of the last equation may be written as 

pik 

! (0); 



uA k = -p-(p + p)Q 



From the Gibbs equation 



Tds = d— + pd— , 

n n 



where s is the entropy per particle, we obtain 



(8) 



(9) 



p = nTs + (p + p) — ■ 
n 



Using here the particle number balance (p) yields 



u i T (0);k = - nTk -(P+P) r 



Inserting the last relation into the balance (0) the latter becomes 



P+P T _ 1 
T 2 



Taking into account the standard definition 



S l = 



rpik rpik 

1 ~ J (o) 



Sl 

T 



;k 



of the entropy flow vector S l of the Eckart theory allows us to write Eq. (Q2h as 



— nsT 



P + P T _ 1 
T 2 

P + P 



rjilk rpik 



r — \T l 

T 2 



-^(0) ) £m/T9ik 



(10) 



(11) 



(12) 



(13) 



(14) 



where £ Ui /T denotes the Lie derivative with respect to the temperature vector u a /T. For T = equation JT4) reduces 
to the well-known expression for the entropy production density of a fluid with constant particle number |7|,p0[|. In 
the latter case one has vanishing entropy production for £ Ui /T9ik = 0, i.e., if u a /T is a KV. This statement has a 
counterpart on the level of kinetic theory following from the requirement that the equilibrium distribution function 
for a simple gas obeys Boltzmann's equation (see, e.g, (_]"§] and below). 

For r 7^ 0, however, there exists the following possibility to characterize a different kind of equilibrium states in 
case u a /T is a CKV, i.e., if 



is fulfilled. Eq. ([lij) then becomes 

S\i - nsT 

Since the CKV property implies the relation 



P + P ] 
T 
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(15) 



(16) 



(17) 



the right-hand side of Eq. (Wm vanishes under the condition 



(P + P) 







(18) 



The entropy production density in such a case is given by 
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&i = nsT , 



(19) 



where T > 0, i.e., production of particles is required to satisfy the second law of thermodynamics. The term nsT 
on the right-hand side of Eq. (|l^) describes the increase of the entropy of the system due to the fact that its phase 
space is enlarged. Each of the created particles carries the same entropy s as the particles already around. From 
the moment of its creation each particle contributes the equilibrium amount s to the entropy of the system. This 
is equivalent to the statement that the produced particles share the equilibrium properties of the already existing 
fluid particles from the very beginning. It is only the amount of entropy describing the very existence of additional 
fluid particles at equilibrium with the particles already there which counts on the right-hand side of Eq. (|l^) . In a 
microscopic picture, all particles including the just produced ones are governed by the same equilibrium distribution 
function (see Eq. (f73|) below). Especially no terms describing dissipative processes such as viscosities or heat fluxes 
appear on the right-hand side of Eq. (|l9|). This justifies the classification of the relationship ( |l9| ) as an equilibrium 
situation although there is an increase in the entropy of the system. We define the "generalized equilibrium " by (i) 
the vanishing of S a a — nsT, i.e., by Eq. with T > together with (ii) the CKV property ( yl| ) of ui/T. For 

r = this requirement reduces to the condition ir — (cf. Eq. (jl§|)), a case dealt with by Bedran and Calvao 
Equation (|l^) for V ^ 0, however, coincides with the frequently used condition for "adiabatic particle production" . 
This condition is a consequence of the definition ( |l3| ) and the generalized equilibrium requirements. There is no need 
to postulate the existence of a supplementary "creation pressure" separately. The existence of such kind of pressure 
is a consequence of the equilibrium conditions and needs not to be assumed. 

To obtain an alternative expression for the entropy production density we realize that differentiating the expression 
©yields 

ai r • , V a q a q a (V a T . \ 

S-i ~ nsT = ns+ —— + — —— + u a , (20) 



T T \ T 

where V a <7 a = h b a q a b etc. On the other hand, the energy balance (Q) may be written as 

p = -6 (p + p + 7r) - V a q a - 2u a q a - o- ab n ab (21) 

with 

(Tab = ^ ^V Q U{, + VbU a - ^h ab Q^j . (22) 

Consequently, via Eqs. (|To| ) and (JsJ) the quantity s is given by 

nTs = -On -{p + p)Y- V a q a - 2ii a q a - a ab 7r ab . (23) 



Combining relations (BOJ) and (p3J) we obtain 



which is an alternative way of writing Eq. ( J14| ) ■ Since the CKV conditions ( |15[ ) imply 



T 



Ua = (Tab = , (25) 



the right-hand side of Eq. (g4|) vanishes, provided condition ( |l8| ) is satisfied. It may be worth mentioning that 
condition (|l^) is not generally equivalent to s — but only if the heat flux q a on the right-hand side of Eq. ( p3| ) does 
not contribute. An important consequence of the CKV property ( |l5| ) is that it completely fixes the evolution of the 
temperature: 

f 6 

- = -- • (26) 
T 3 K ' 

On the other hand, the behaviour of the temperature is determined by general thermodynamical considerations. The 
fluid equations of state may be written as 

P = p(n,T) , P = P (n,T) . (27) 
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Differentiating the latter relation and using the balances (Bh and ( |lC)| ) we obtain 



dp _ (Sp/dT)„ dp _ (dp 



where the abbreviations 

<9p ~ {dp/dT) n ' <9T _ V<9T, 

have been used. Comparing Eqs. ( p6|) and (^8| ) and assuming fom now on q a = 0, i.e. s — 0, yields 



3 dp 



9 . (29) 



The CKV property fixes the creation rate. An equilibrium according to S? n — nsF = (generalized equilibrium) is 
only possible if fluid particles are created at a specific rate given by Eq. (E9). The production rate F is positive for 
any equation of state < p < p/3. It vanishes for p = p/3, i.e. for radiation. This means recovering the well-known 
fact that a fluid obeying the equation of state p = p/3 may be in "global equilibrium" if the spacetime admits a 
timelike CKV. The CKV property of Ui/T implies the temperature law (^6|) which just coincides with the well-known 
temperature behaviour (adiabatic cooling) of a relativistic fluid in the expanding universe. For r = only a fluid with 
the equation of state p — p/3 satisfies Eq. (p6|). For any other equation of state the equilibrium conditions require 
(f) = in Eq. (|l5|), i.e., the existence of a KV. If one drops the condition F — 0, i.e., if one allows the particle number 
to increase, radiation is no longer a singular case. Generalized equilibrium states (|l5|) with <fi ^ are possible for any 
equation of state in the range < p < p/3. For any of these equation of state there exists a specific production rate 
( p9| ) which is necessary to satisfy the generalized equilibrium conditions (]l5|) and ([0]). The production of particles 
at a rate different from Eq. (E9T) disturbs the generalized equilibrium. In other words, a fluid with an equation of 
state in the range < p < p/3 may only be in generalized equilibrium if fluid particles are produced at a specific rate 
depending on the equation of state. A specific rate F on the other hand, uniquely fixes tt acording to condition jl^): 

* = -(P + P)(l-l¥) ■ (30) 



3 dp 

We recall that the generalized equilibrium conditions naturally imply the existence of an effective bulk pressure, i.e., 
the latter has not to be postulated separately. It is obvious that tt as well as F vanish for p = p/3, i.e. for radiation. 
For nonrelativistic matter, characterized by equations of state p = nT and p = nm + \nT with m^> T one obtains 

r = ie, n^-^p + p)^-^ , (m»T), (3i) 

i.e., r is half the expansion rate and \ir\ is a remarkable fraction of the energy density. 



III. THERMODYNAMICS OF THE CREATION PROCESS 



With equations of state of the general form (27) the particle number density n and the temperature T are the 
primary thermodynamic variables of our system. While the (generalized) equilibrium behaviour of the temperature 
is given by Eq. ( |26| ) independently of the equation of state, the particle number density changes according to 



n 
n 



ldp 
3 dp 



e , 



where we have combined Eqs. (|) and Use of Eq. @) in Eq. @ yields (recall that q a = a ab = 0) 



P= -^ {P + P) dp- 



(32) 



(33) 



Differentiating the first equation of state ( |27| ) and applying relations (J6|) and (|28|) under the condition s = 0, i.e., with 
the equations ([l8]) and (p9|), one finds 
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p = c s p 



(34) 



where 



dp 

d P J ad 



n dp T (dp/dT) 2 



p dn 



p dp/dT 



(35) 



is the square of the adiabatic sound velocity. An essential quantity to characterize the particle production process is 
H/T, where fi is the chemical potential. Using Eqs. (pa) and (p4[) in the Gibbs-Duhem relation (see, e.g., 0) 



dT 



d P =(p + p) — + nTd[^ 



provides us with 





e p + p 


dp_ 




p _ dp' 


e 




3 nT 




p dp 


"3 



(36) 



(37) 



Here we have used that for p = nT, i.e. for a classical gas, the square of the sound velocity ( pq ) may be written as 

9 nT 



p + p 



dp 



(38) 



The right-hand side of Eq. (|37j ) vanishes only for radiation (c 2 s = 1/3 and dp/dp = 3). Characterizing the equilibrium 
states under consideration also requires information about the spatial derivative of pt/T . From the momentum balance 



(p + p + TT) u m + V m (p + 7r) = 



(where we assumed q a — n — 0) and the Gibbs-Duhem equation (p6J) together with Eq. (|25|) we obtain 



-1TU a 



V Q 7T 



(39) 



(40) 



where n is given by Eq. (|30[). It follows that also the spatial derivative of p/T vanishes for radiation. Restricting 
ourselves to a classical gas and using Eqs. (|(j|) and (P^) , Eq. (^) may be written as 



VqTT 



=- l-r^ 



1 3p\ V a T nT 



3 dp J T 



P- 



(41) 



On the other hand, an alternative expresssion for V a 7r in terms of V a T and V a /x may be obtained directly from Eq. 
©. To this purpose we first realize that the Gibbs-Duhem equation ([56]) implies 



nT 
p + p 



With p = nT we obtain 



From the Gibbs equation @ it follows that 



nT V a n 
-V a s H 



p + p p + p 

Spatial differentiation of the second equation of state (p7j) provides us with 



VaP 

P + P 

where we have used the equation 



1 



T dp\ V a n T 9pV a r 
p + p J 71 P + P dT T 



(42) 



(43) 



(44) 



(45) 



6 



dp p + p T dp 
dn n n dT ' 



(46) 



which follows from the fact that the entropy is a state function. Comparing the expressions (|44|) and ([15|) and using 
Eq. (H) yields 



V Q is 



Realizing now that || 

use of p — nT and Eqs. ([l3|) and 

V Q 7T 



TJ ~ [dp p\ T 



u _ p + p 
S+ T nT ' 

yields for the spatial gradient of the effective viscous pressure n in Eq. 



(47) 



(48) 



p + p 



= - 1- 



ldp 
3 dp 



- l- 



id P \ p + pV a T 



3 dp J nT T 



nT 
p + p 



ldp\ (dp p\ V a T 1 



3 dp J \dp p ) T 



dp 



For a gas one has ul 



dp ( m \ 
~dp~ ~ \ f~) 



2 -l + 5 P±P 
nT 



P + P 
nT 



(49) 



(50) 



Consequently, the expression ( (49| ) for V a 7r becomes, in terms of V a T and V a (p/T), 



P+P 



I dp 



3 dp 
nT 
p + p 



= 1-77 V a £ - 1--^ 



ldp\ p + pV a T 



TJ V SdpJ nT T 
1 dp\ f dp p\ 2 /m\ 2 1 



, p + p \ f^p_p 

nT ) \ dp p 



T 



(51) 



Comparing now the expressions (|4l|) and (pl|) provides us with a relation for V a (p/T) in terms of V a T: 



-^)V.(£] = 



p + P 3 dp 



T 



nT 



3 dp ) \p p + p 



dp_ _ p 
dp p 



2 /m\ 2 

3 \T> 



2p 



3 P ) \dp p 



dp p 



T 



(52) 



The case of vanishing parenthesis in front of V a (p/T) on the left-hand side of the last equation, i.e., pj {p + p) — 
(1/3) dp/ dp — corresponds to (cf. Eq. (30)) tt = — p = —nT. In such a case the momentum balance ( p9| ) reduces to 
u m = which, according to Eq. (pa), implies V a T = 0. Then, V a (p/T) by Eq. (p3|) is simply given through 



'•( 



(53) 



Of special interest are again the limiting cases of radiation and nonrelativistic matter. In the first case, i.e., for 
p = p/3, Eq. © yields 



= 



(p = p/3) 



(54) 



The quantity p/T is constant both in space and time (cf. Eq. (|37|)). In the opposite limit p — nT and p = nm + \nT 
with T one obtains 



V 



q (t) 



m VqT 
T T 



(m> T) 



(55) 



which may be written as 
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V Q 0^-) =0 (m » T) . (56) 

Since one finds from relation (|37j ) that also 

{^f) -0 (m»T) (57) 

holds in this limit, it follows that the quantity (p, — m) /T is generally constant for nonrelativistic matter. It is 
well-known (see, e.g., |5|,|2^]) that // — m is just the nonrelativistic chemical potential. 

This completes our general discussion of the thermodynamics of fluid particle creation in generalized equilibrium. 
We are now prepared to discuss the implications of this kind of particle production for the cosmological dynamics. 
While part of this discussion recalls earlier work one should keep in mind the much more general level of the present 
considerations. While the discussion in jl2| relied on a very specific model for the creation process, the conclusions 
of this paper are completely model independent. 



IV. PARTICLE PRODUCTION AND COSMOLOGICAL DYNAMICS 



The main advantage of a fluid approach to particle production is the possibility to calculate the backreaction of 
this process on the cosmological dynamics. As was shown in sections II and III this problem is equivalent to studying 
the dynamics of a bulk viscous fluid universe. Moreover, the effective bulk pressure that describes the production 
process in the equilibrium case is explicitly given in terms of p and p. implying the dynamical laws (p2h and (pq). 
The expression dp/ dp in the latter equations is generally given by Eq. (|50|). Explicit integration of ( p2[ ) and (|33|) 
is possible in the limiting cases of radiation p = nT, p = 3nT and nonrelativistic matter p — nT, p = nm + |?iT, 
m 3> T. Introducing a length scale a according to 

6 = 3^ (58) 

one recovers the familiar dependences 

n oc a~ 3 , p cx a~ 4 (p — p/3) (59) 

for radiation since T as well as 7r vanish in this case. For nonrelativistic matter on the other hand, one obtains 

n cx a~ 3/2 , p cx a _3/2 (m > T) . (60) 

The temperature behaves according to Eq. ( |26| ) for any equation of state. Since nonrelativistic matter is produced at 
a rate T = 0/2 (cf. Eq. (|l])) its standard behaviour for the case of conserved particle number ( n cx a~ 3 , T cx a~ 2 , 
p cx a~ 3 ) changes considerably. All the thermodynamic quantities decrease more slowly now since the decay of n, T 
and p due to the expansion is counteracted by corresponding production terms. The differences in the behaviour of 
n, T and p acoording to the relations ( |60| ) compared with the standard behaviour for T = are consequences of the 
backreaction of the production process on the fluid dynamics. Also the expansion of the universe is modified in case 
matter in generalized equilibrium becomes dynamically dominating. In a homogeneous and isotropic universe the 
lenght scale a coincides with the scale factor of the Robertson- Walker metric. Restricting ourselves to the spatially 
flat case, the scale factor obeys the equation 

a 2 

3— = Kp , (61) 
a z 

where k is Einstein's gravitational constant. For radiation we recover, of course, a cx t 1 / 2 . Inserting, however, the 
energy density p from Eq. (|6(]) into Eq. ( |6l| ) we find that the scale factor behaves such as 

a cx t 4/3 (m»T) , (62) 

i.e., a > instead of the familiar a cx t 2 / 3 with a < for p cx a~ 3 corresponding to T = 0. The production of massive 
particles in equilibrium implies accelerated expansion ( a > 0), i.e., power law inflation. In other words, only in a 
power law inflationary universe an equilibrium such as discussed in this paper is possible. 
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The radiation and nonrelativistic matter equations of state constitute the limiting cases of the fluid behavior. The 
backreaction is largest for m^> T and it vanishes for m <C T, i.e., m — > 0. There exists an intermediate equation of 
state yielding a cut with a — 0. The condition for o = is 

p + 3 (p + tt) = , (63) 

or, with 7r from Eq. ((30|), 

dp 2/3 



<9p p + p 

Inserting here Eq. (p0[), we obtain the following cubic equation for h/T = (p + p) / (nT): 



(64) 



h\ 3 (h^ 



T) 5 \f 

The quantity h is the enthalpy per particle, given by 



vf) 



--2 = 0. (65) 



7>^ /m\ 

ft = m ±l^xl . ( 66 ) 

K% and ^3 are Bessel functions of the second kind. The critical value (m/T) cr that corresponds to a = turns out 
to be (m/T) cr « 9.55. We have a > for (m/T) > (m/T) cr , while (m/T) < (m/T) cj , corresponds to a < 0. In 
terms of a "7" law, i.e., p = (7 — 1) p, the corresponding critical value "f cr is 7 cr ~ 1.09. There is power law inflation 
for 1 < 7 < 7 cr . The case 7 = j C r corresponds to a = 0. We mention again that all the results concernig the 
backreaction of equilibrium particle production on the cosmological dynamics are quite general and do not rely on a 
specific creation model. 

In the following section we will establish relations between the results obtained so far and corresponding con- 
siderations on the level of relativistic kinetic theory. This will provide us with a microscopic justification of the 
phenomcnological generalized equilibrium concept. 



V. KINETIC THEORY AND PARTICLE PRODUCTION 



A. General relations 



For conserved particle numbers the fluid dynamics of a gas may be derived from Boltzmann's equation for the 
one-particle distribution function. One may address the question whether a corresponding microscopic justification 
also exists in the case of varying particle numbers. Previously |n],[^j a modified Boltzmann equation was proposed 
in which an additional "source term" describes the change of the one-particle distribution function due to particle 
number nonconserving processes, supposedly of quantum origin. 

The one-particle distribution function / = / (x,p) of a relativistic gas with varying particle number is supposed to 
obey the equation 

L [/] = p'f,i -r^V|4 =C[f]+H (x, P ) , (67) 

where / (x,p) p k nkdT,dP is the number of particles whose world lines intersect the hypersurface element n^dS around 
x, having four- momenta in the range (p,p + dp). 

dP = A(p)S (p l Pi + m 2 ) dP4 is the volume element on the mass shell p % Pi = —m 2 in the momentum space. A(p) = 2, 
if p l is future directed and A(p) = otherwise; dP^ = — gdp° dp 1 dp 2 dp 3 . 

C[f] is the Boltzmann collision term. Its specific structure discussed e.g. by Ehlers [^) will not be relevant for our 
considerations. Following Israel and Stewart || we shall only require that (i) C is a local function of the distribution 
function, i.e., independent of derivatives of /, (ii) C is consistent with conservation of four-momentum and number 
of particles, and (iii) C yields a nonnegative expression for the entropy production and does not vanish unless / has 
the form of a local equilibrium distribution (see Eq . ( |73| ) below) . 

The term H(x,p) on the right-hand side of ( |67|) takes into account the fact that the number of particles whose 
world lines intersect a given hypersurface element within a certain range of momenta may additionally change due to 
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creation or decay processes. On the level of classical kinetic theory we shall regard this term as a given input quantity. 
Below we shall give an example for the possible functional structure of H(x,p). 

By the splitting of the right-hand side of Eq. ( |67| ) into C and H(x,p) we have separated the collisional from the 
creation (decay) events. In this setting collisions are not accompanied by creation or annihilation processes. In other 
words, once created, the interactions between the particles are both energy-momentum and number preserving. For 
a vanishing H(x,p) Eq. ( |67|) reduces to the familiar Boltzmann equation (see, e.g., [[|-f|||]). 

The particle number flow four-vector N l and the energy momentum tensor T zk are defined in a standard way (see, 
e.g., §) as 



n 1 = I dPff( x , P ) , r fe = / dPtfpKf(x, P ) 



(68) 



While it will turn out that in the equilibrium case the first moment of / in Eq. (pq ) may be identified with the 
quantity (|4|), we have used the symbol N l immediately. The second moment of / denoted by T ik in E q. ( | 68|) will not, 
however, coincide with the energy momentum tensor T lk in Eq. ([!]) with Eq. (^|). The integrals in (|68f) and in the 
following are integrals over the entire mass shell p l pi = —to 2 . The entropy flow vector S a is given by 




S a = -Jp a [/In/-/] dP, 



(69) 



where we have restricted ourselves to the case of classical Maxwell-Boltzmann particles. 
Using the general relationship [Q 



p ai ....p a "p b fdP 



= J p a \..p a "L[f]dP 



and Eq. (67) we find 



K= (C[f]+H)dP, f ak = p a (C[f} + H)dP , 



and 



S a n 



J ln/(C[/] 



H) dP 



(70) 



(71) 



(72) 



In collisional equilibrium which we shall assume from now on, In / in Eq. ([72|) is a linear combination of the collision 
invariants 1 and p a . The corresponding equilibrium distribution function becomes (see, e.g., Q|) 



(73) 



(74) 



f°(x,p) = exv[a + [3 a p a ] , 

where a = a (x) and j3 a (x) is timelike. Inserting the function ( |73| ) into Eq. ([37]) one obtains 

[p a a, a + (3 {a ., b)P a p b ] f° = H(x,p) . 

It is well known |&0 that for H(x,p) ~ this equation, which characterizes the "global equilibrium", admits 
solutions only for very special cases in which a = const and /3 a is a timelike Killing- vector. The only exception is the 
case m = 0, in which j3 a is only required to be a CKV in order to satisfy the equilibrium conditions. It will be our 
main objective in this section to explore whether a nonvanishing source term H (x,p) admits Eq. ( [74|) to be fulfilled 
under the condition that f} a is a CKV also for m > 0. With the equilibrium distribution funktion (|T3|) the balances 
Til) reduce to 



N? a = J HdP , f ak k = J p a HdP . (75) 
In collisional equilibrium there is entropy production only due to the source term H . From Eq. ( f72|) we obtain 

(76) 



s -a = -J H{x,p) in Alp 



implying Sf a = -aN« a - (3 a T ab b . With / replaced by f° in Eqs. (f|) and (|g), N a , T ab and S a may be split with 
respect to the unique four-velocity u a according to 
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N a = nu a , T ab = pu a u b +ph ab , S a = nsu a . (77) 

The exact integral expressions for n, p, p and s are given by the formulae (177) - (180) in g. 
Using the first Eq. fin\ ) and defining 

T = - [ H(x,p)dP , (78) 
n J 

the first Eq. (|75|) coincides with Eq. (^). It is obvious that T is the particle production rate. Similarly, with the 
decomposition (|77|) and the abbreviation 

t a = - J p a H(x,p)dP , (79) 

the second equation (f75|) may be written as 

f.f + t a = . (80) 

For t a ^ the energy momentum tensor T ak is not conserved. We are interested in the question whether it is possible 
to replace the left-hand side of the last equation, T a k k + t a , by a locally conserved quantity, T ak , 1 in such a way that 
the source term t a in Eq. (|80| ) is mapped onto an effective viscous pressure ir of the energy momentum tensor T ak , 
which is to be identified with the corresponding quantity in Eq. (Q). 
The equations ( |so|) imply 

p + <d(p + p) =u a t a (81) 

and 

(p + p) u a + V aP = -h ai t l . (82) 
Comparing these balances with relations ( ^l| ) (for q a = a a b = 0) and ( |39"| ) suggests the tentative identifications 

9tt = -u a t a (83) 

as well as 

nu a + V a 7r = h ai f . (84) 

In order to check the consistency of this interpretation it is necessary to specify the source term H(x,p) which 
describes the variation of the one-particle distribution function due to the change in the number of gas particles. In 
the following subsection we will generalize previously used expressions for H (x,p). 

B. Effective rate approximation 

The quantity H (x,p) is an input quantity on the level of classical kinetic theory. It is supposed to represent the 
net effect of certain quantum processes with variable particle numbers (see, e.g., [|l8|,|l9|). Also internal interactions 
within the cosmological fluid may give rise to an increase in the overall number of particles [ p^|Jl^ ] . Following earlier 
considerations []ll"| , [l2"[ we assume H (x,p) to depend linearly on the equilibrium distribution function (|73|), i.e., 

H(x,p) = C(x, P )f°(x,p) . (85) 

We also assume that £ (x,p) depends only linearly on the momenta p a . This ensures that the source terms T and t a 
in the balances (^) and ( |80| ) do not depend on moments of the distribution function higher than the second one. The 
four-momenta p a may be decomposed into 

p a = Eu a + Xe a , (86) 

where e a is a unit spatial vector, i.e., e a e a = 1, e a u a — 0. Consequently, one has E = —u a p a and A = e a p a and 
the mass shell condition p a p a = —m 2 is equivalent to A 2 = E 2 — m 2 . Moreover, /i a bP a p b = A 2 is valid. These 
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properties suggest assuming £ to consist additively of three terms: A term v which does not depend on p a at all, a 
term proportional to E = — u a p a , and, finally, a term proportional to A = e a p a - 



C(x,p) = v{x) + 



A 



E 



t (x) a (x) 



(87) 



This ansatz for Q represents the most general structure compatible with a linear dependence of £ onp°. The rate of 
change of the (equilibrium) distribution function is then determined by the three spacetime functions r, and a. 
Different from earlier considerations pd| , ^2| we have introduced here a term proportional to A = e a p a = \J h a bp a p b . 
As we shall see below this more general structure will remove some unpleasant properties of the model studied in 

OS 

With the relations (j8q) and (p7|) for H (x,p) the source terms nT and t a in Eqs. (f78[) and (J79|) are given by 



nT = v (x) M (x) - —N a — vM + — , 

T T 

where M is the zeroth moment of the distribution function, M = J dPf (x,p), and 

u b T ab e b T ab 



t a 



respectively. The last relations imply 



-vN a 



. P 
vn H — 

T 



t (x) a (x) ' 



h C it l = p ■ 



(88) 



(89) 



(90) 



Via the functions is, r, and a the source terms T, Uit % , and h C it l are coupled to the fluid quantities n, p, and 
p, representing the first and second moments of the distribution function. Because of the restriction to a linear 
dependence of £ on p a in Eq. (87), no higher-order moments appear in the expressions ( |38| ) and (p9|). Using the 
structure (Efj) with the ansatz (87) for the source term H(x,p) in Eq. (|7^), the latter reduces to the generalized 
equilibrium condition 



With the decomposition 



P a a,a + f3{a;b)P a p h = 

condition ( |9l| ) is satisfied for 



E 

T 



A 
a 



1 

a = — 

T 



e Q V a a = - 
a 



and 



with 



P(a;b) = H x )9ab , 

y/m 2 . Together with Eq. ( |l7j ) this implies 



ie 

3T 



(91) 



(92) 



(93) 



(94) 



Using the identification j3 a = Ug/T we have recovered the CKV property ( [l5| ) as a (generalized) equilibrium condition 
for the distribution function ([73|). According to the relations (|9^) the function a may vary both in time and in space. 
The possibility of a spatial variation is due to the last term in the expression (87) and was not taken into account in 
As we will see below, a spatially varying a is, however, essential to establish the full correspondence between 
the microscopic and macroscopic levels of description. Identifying a with p/T of the phenomenological theory, the 
relation (76) may be written as 



na 



P r, U a t a 
Til — 



T 



T 



Using here Eq. (Eq) we obtain 



S a „ - nsT = - 



P + P , 
T 



Ugt a 

T 



(95) 



(96) 
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Introducing the abbreviation (|83j) for ir into the last relation and taking into account the conditions ( |25| ) we recover 
Eq. (H). For q a = the right-hand side of Eq. © coincides with nTs (cf. Eq. @). With a = fi/T, the first 
condition (B2h may be combined with Eq. (B7n to yield 



r 3 dp, 



(97) 



On the other hand, r 1 and T are related by Eq. (p8[). It follows with Eq. (|29| ) that the function v is given by 



JL(i_£.|e. 



The second equation (193) together with Eq. (43) fixes the function a 



(98) 



v(£)v.(£) 



V a n pV a T 



p r 



n p T 



(99) 



An alternative expression for 1/a maybe obtained by using formula ( p2| ) for V a (p/T). 

With the first relation (^) and the expressions ( |9^ ) for and (|97j) for r _1 the energy balance (|8l"l ) is consistent with 
the balance p + ® (p + p + ir) = of the phenomenological approach with ir given by relation (|30|) , equivalent to Eq. 
(33). To check the corresponding property for the momentum balance one has to show the equivalence between Eqs. 
(32) and j39|). The identification a = p/T allows us to combine Eq. ( jif^ ) with the second equation j92]). Together 
with the second relation (Q) this justifies relation (|84|). Consequently, the balances ( p0| ) with a nonconserved energy- 
momentum tensor T ak given by T ab = pu a u b + ph ab are entirely equivalent to the local conservation laws T a \ = for 
the energy-momentum tensor T ak — pu a u k + (p + tt) h ak . The source terms t a in the balances ( |S0| ) have consistently 
been mapped onto the quantity it in the conserved energy- momentum tensor T ak . This completes our microscopic 
derivation of generalized equilibrium including the equivalence between particle production processes and effective 
viscous pressures. 



VI. SUMMARY 



The "generalized" equilibrium concept introduced in this paper allows one to characterize cosmological fluids with 
growing number of particles. This equilibrium relies on the conformal Killing vector property of Ui/T and implies 
nonvanishing entropy production which, however, is entirely due to an enlargement of the phase space of the system 
but not to dissipative processes. For a gas the "generalized" equilibrium may be derived from relativistic kinetic 
theory The corresponding gas particles are governed by an equilibrium distribution function although their number 
is generally increasing. In the limit of relativistic particles with equation of state p = p/3 the "generalized" equilibrium 
reduces to the well-known "global" equilibrium of standard relativistic thermodynamics. Nonrelativistic particles may 
be in generalized equilibrium only in a power law inflationary universe. The latter statement confirms and generalizes 
earlier results. It was shown here to be model independent. 
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